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Abstract 

In this paper, a method to construct topological template in terms 
of symbolic dynamics for the diamagnetic Kepler problem is proposed. 
To confirm the topological template, rotation numbers of invariant 
manifolds around unstable periodic orbits in a phase space are taken 
as an object of comparison. The rotation numbers are determined 
from the definition and connected with symbolic sequences encoding 
the periodic orbits in a reduced Poincare section. Only symbolic codes 
with inverse ordering in the forward mapping can contribute to the 
rotation of invariant manifolds around the periodic orbits. By using 
symbolic ordering, the reduced Poincare section is constricted along 
stable manifolds and a topological template, which preserves the or- 
dering of forward sequences and can be used to extract the rotation 
numbers, is established. The rotation numbers computed from the 
topological template are the same as those computed from their orig- 
inal definition. 
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1 Introduction 



Many interesting nonlinear systems in experiments are wells described by 
low-dimensional dynamical models. Their dynamical processes include of 
transition of periodic orbits from stable to unstable and bifurcation to chaos. 
To make a global understanding of the system, some topological and geomet- 
rical methods based on periodic orbits are developed]!]. Symbolic dynamics, 
as a coarse-grained description of the dynamics, provides an effective tool to 
depict the topological dynamics[2]. In special, for a chaotic system, two curve 
families of stable and unstable manifolds intersect each other and decompose 
a Poincare section[3,4]. Enumeration and existence of unstable periodic or- 
bits (UPOs) in the Poincare section are determined [5, 6] and some numerical 
methods such as finding UPOs are proposed [7]. 

Besides the features of dynamics in a Poincare section, evolution of man- 
ifolds in a phase space, which possibly gains a geometric insight into the 
dynamics, is another important characteristic of the system. For example, 
bifurcation to chaos can be identified by tracing the evolution of stable and 
unstable manifolds in the phase space with parameters. Recently, some al- 
gorithms for computing two-dimensional stable and unstable manifolds in 
a three-dimensional phase space are proposed[8,9,10] and applied to visual- 
ize the structure of chaos[ll]. Since UPOs are closely related to invariant 
manifolds, a basic issue in their relations is describing how local invariant 
manifolds rotate around UPOs in the phase space. Moreover, how reducing 
their topological relation into a two-dimensional template is important for 
understanding the global organization of UPOs in chaos. 

To motivate the visualization of rotation of invariant manifolds around 
UPOs and construction of topological template in reduction of stable man- 
ifolds, we consider the model for diamagnetic Kepler problem (DKP)[12]. 
In our previous works, two coordinate axes are chosen as a Poincare sec- 
tion to form an annulus in a lifted space. The dynamics on the annulus 
can be reduced by considering the symmetry of system. In view of stretch- 
ing and wrapping in the lifted space, symbolic dynamics without involving 
bounces has been established [13]. Due to the ordering of stable and unstable 
manifolds in the minimal domain (a reduced Poincare section) , a method to 
extract UPOs corresponding to short symbolic strings is proposed. A one to 
one correspondence between UPOs and symbolic sequences is shown under 
the system symmetry decomposition [14]. Although we only focused on the 
case of zero scaled energy, in our numerical experiences, the methods can be 
still used at the scaled energy ±0.1. In this paper, for the DKP, we calculate 
rotation numbers of invariant manifolds around UPOs and set up a con- 
nection between rotation numbers and symbolic sequences. Using symbolic 
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dynamics, we reduce stable manifolds to construct a topological template, 
which preserves the topological relation of UPOs. 

2 Model system and Poincare section 

The Hamiltonian of a hydrogenic electron (with zero angular momentum) in 
a uniform magnetic field B directed along the z-axis is given by 

H = (l/2m)(p 2 p + p\) - e 2 /(p 2 + ^ 2 ) 1/2 + \nu? p\ (1) 

where u = eB/2mc is half the cyclotron frequency. Converting to atomic 
units and transforming the cylindrical coordinates to semiparabolic ones, the 
Hamiltonian becomes 

h = f + f - e(/x 2 + v 2 ) + ^AV + v 2 ) = 2, (2) 

where e = E^~ 2 ^ is the scaled energy depending on the energy E and the 
dimensionless field strength parameter 7 = 2ou. The symmetry group of 
h consists of the identity e, two reflections a^, a u across the /i, v axes, two 
diagonal reflections o"i 3 , <r 2 4, and three rotations C 4 , C 2 and Cf by ir/2, n and 
37r/2 around the center, respectively [15]. In the following, the symmetries 
C4, C2 and Cf are denoted by p, n and p, respectively. The time- reversal 
symmetry is denoted by T. 

Figure 1 displays an orbit and boundary of the transformed potential for 
e = 0. A Poincare section is chosen as follows. Imagine that the p and v 
axes are both of a finite width and length. A counter-clockwise contour is 
taken along the perimeter of the area forming by the two crossing imaginary 
rectangles. The Poincare section is then obtained by recording the position 
and the tangent component of the momentum along the contour, i.e., the 
Birkhoff canonical coordinates[16] at intersecting points with the contour 
where an orbit enter the inside of the contour. The length of the contour 
is infinite. It is more convenient to transform the contour to one with a 
finite length. For example, in the first quadrant, the transformations s = 
— along the positive \i axis and s = v/(l+v>) along the positive v axis 
convert the segment of the original contour in the first quadrant to interval 
of length 2 parametrized with s G [—1, 1). The variable corresponding to the 
momentum is taken as v — —p^/p at the positive \i axis, and v = p v /p at the 
positive v axis, where p = ^p 2 ^ + p*. In this way we may parametrize the 
whole contour with s G [—1,7) and define corresponding v. The rotational 
symmetry under p, n and p in the original configurational space becomes the 
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translational symmetry of shifting s by a multiple of 2 in the s — v plane. 
The dynamics on the Poincare surface is then represented by a map on the 
annulus s G [—1,7) and t> G [—1,1], which is taken as a fundamental domain 
(FD). When consider an image and preimage of the FD, we need extend to 
its lifted space. The partial image and preimage of the FD in the lifted space 
is given in Fig. 2. For example, in Fig. 2(a), zones 1 (2) and 1' (2') in the 
strip 1 (2) are mapped forward into zones +3 (+0) and +2 (+1), respectively. 
In the same way, the backward mapping of the strip 1 is given in Fig. 2(b). 
Since the rotational symmetry of the Hamiltonian (2) corresponds to the 
translational symmetry in lifted space, the annulus (s G [—1, 7), v G [—1, 1]) 
on the Poincare section can be reduced to a domain (s G [0, 2), v G [—1, 1]). 

In the conservative system, classical dynamics preserves an invariant vol- 
ume in the phase space under constricting, stretching and folding. This 
behavior can be displayed in the Poincare section. In Fig. 3(a), we draw the 
9 lines (s G (0, 1), v G [—1, 1]) in the reduced domain (RD) and their forward 
mapping in the lifted space. In the mapping, the original zone in the RD 
is stretched and folded, as well as wrapped. In order to display ordering of 
the lines in the forward mapping, we also plot connecting lines between two 
different strips in the lifted space. In the lifted space (s G (2, 3), v G [—1, 1]), 
the ordering of lines in the top-left part preserves the same as the original 
one and the ordering in the bottom-right part is in reverse. In another lifted 
space (s G (1,2), v G [—1,1]), the ordering of lines in the whole part is 
in reverse. It is clear that the wrapping of lines in the forward mapping is 
clockwise, if the lines in the strip (s G (1, 2), v G [— 1, 1]) are stuck to those in 
the strip (s G (2, 3), v G [—1, 1]) in terms of their ordering. In the same way, 
in Fig. 3(b), the similar result can be obtained from the 9 lines (s G (1,2), 
v G [—1, 1]) in the RD and their forward mapping in the lifted space. So, the 
forward mapping illustrates the rotation of RD in the clockwise direction. 

3 Stable and unstable invariant manifolds 

In general, the invariant manifolds as a subset are contained in manifolds. 
The method for calculating stable and unstable manifolds (dynamical folia- 
tions) used in through the rest of our works is detailed in [2] . Using the same 
method, we can thus generate the stable and unstable invariant manifolds 
through unstable periodic points in a two-dimensional Poincare section. Here 
we present a short introduction of the method. 

(i) Unstable manifolds: Taking a circle around n steps backward mapping 
(x- n , D-n) of an unstable periodic point (xo, yo), we get an ellipse centered 
at (x , y ) after the same steps forward mapping. Its long axis points to the 
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most stretching direction. When we fix the point (x , y ) and increase n, 
the ellipse is stretched and rotated, as well as its most stretching direction 
changes slightly. When n — > oo, the most stretching direction approaches a 
limit. This direction is the most stretching direction of the point (xq, yo)- 
After going a short distance along the direction, we get a new point (xo, yo). 
Repeating the above process, we can get the new most stretching direction. 
Finally, an unstable invariant manifold is generated by connecting the points 
(xq, yo)- In fact, the deformation of ellipse is closely related to the dynamical 
matrix of Poincare mapping. 

(ii) Stable manifolds: Taking a circle around n steps forward mapping 
(%n, Un) °f an unstable periodic point (x , y ), we get an ellipse centered at 
(x , yo) after the same steps backward mapping. Following the above similar 
process, we can get the most stable direction of the point (x , yo) and a 
stable invariant manifold. 

4 Evolution of unstable manifolds in rotated 
Poincare sections 

In order to display the evolution of invariant manifolds around an UPO, 
we rotate counter-clockwise the (/z, v) coordinates to the (/z^, u^) coordi- 
nates with an angle G [0, 7r/2]. By using the same transformations in 
Section II, the Poincare map (s^, v^) is obtained from the (/i^, v^) coor- 
dinates and reduced to a domain (s^ G [0,2), G [—1,1]). Since h has 
Ci v and time-reversal symmetries, we take 4 UPOs (4) (5) (14) (15) with dif- 
ferent symmetries in the Table I as examples to investigate the evolution of 
invariant manifolds. Their plots in the configuration space with 10 rotation 
coordinate axes and periodic points with unstable invariant manifolds in the 
RDs are drown in Figs. 4(a)-(d), respectively. From the initial point in the 
+V4>=o coordinate axis corresponding to the point 1 in the RD with = 0, 
each UPO goes into the second quadrant as displayed in the configuration 
space by an arrow. Its time process is recorded in the RDs with 10 rotation 
angles 0. So, the figure in the RD with = n/2 is the same as that in the 
RD with = besides the first point in the former RD is the second point 
in the later one. In a periodic process, we will calculate the advanced phase 
9 of an unstable manifold in rotation around the periodic orbits to determine 
rotation number 9/2n. The phase 9 is counted positive (negative) when the 
rotation of unstable manifold around the periodic orbits is counter-clockwise 
(clockwise). 

The UPOs with 4 different symmetries are described as follows: 
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(i) The UPO (4) displayed in the configuration space of Fig. 4(a) has 
and T symmetries, but not p, n, p, cr u , cr 13 and (X 2 4 symmetries. In the 
second quadrant of configuration space, the orbit starting from the +z^=o 
coordinate axis goes to the —p<p=o (+^=71-/2) coordinate axis. It corresponds 
to that the point 1 moves to 2 in the RD with = 0, i.e. 1 in the RD with 
= 7r/2. In the process, the orbit passes through the +u ( f, =n /i8, +v<f, =n /g, ■ ■ •, 
+^=47r/9 coordinate axes, while the point 1 moves in the RDs with = n/18, 
= 7r/9, • • •, = An/9. At the same time, an unstable manifold passing 
through the point 1 evolves in the RDs. The phase of unstable manifold 
advances by an angle close to —n. In the third quadrant of configuration 
space, the orbit starting from the — p<f>=o coordinate axis goes to the —v<f,=o 
(—^0=^/2) coordinate axis. It corresponds to that the point 2 moves to 3 in 
the RD with = 0, i.e. 2 in the RD with = ir/2. In the process, the 
orbit passes through the —p<t>=n/i8, —^<t>=n/9, " •, —p<t>=4,n/9 coordinate axes, 
while the point 2 moves in the RDs with = vr/18, = n/9, • • •, = 47r/9. 
The phase of unstable manifold advances by an angle close to — n. The orbit 
starting from the —v<j,=o coordinate axis goes into the fourth quadrant of 
configuration space and back to the — z/0 =o coordinate axis. It corresponds 
to that the point 3 moves to 4 in the RD with = 0. In the process, the 
orbit passes through the — i^ =7r /i8 coordinate axis two times, while the point 
3 goes into the RD with = n/18, then leaps to the point 4 and comes back 
the RD with = 0. In the moving processes of the point 3 from the RD 
with = into the RD with = n/18 and of the point 4 from the RD 
with = n/18 into the RD with = 0, the directions of unstable manifolds 
are almost invariant. However, in the leaping process of the orbit from the 
point 3 to 4 in the RD with = n/18, an unstable direction indicated as 
the arrow from the periodic point to its neighboring point on the unstable 
manifold is approximately reversed. The phase of unstable manifold advances 
by an angle close to —n. In the third quadrant of configuration space, the 
orbit starting from the — z^,=o {—p<f>=n/2) coordinate axis goes to the —p^o 
coordinate axis. It corresponds to that the point 4 moves to 5 in the RD with 
= 0, i.e., the point 3 in the RD with = n/2 moves to 5 in the RD with 
= 0. In the process, the orbit passes through the —^0=4^/9, —p<f>=7 n /i8, 
■ ■ ■, —p<f>= n /i8 coordinate axes, while the point 3 moves in the RDs with 
= 47r/9, = 7n/18, ■ ■ -, = n/9 and then the point 5 moves in the RD 
with = n/18. Since the —v^= n /is coordinate axis intersects the orbit in the 
fourth quadrant, two points 3 and 4 are added in the RD with = n/18. 
The phase of unstable manifold advances by an angle close to — n. In the 
second quadrant of configuration space, the orbit starting from the —p<f>=o 
(+^=,1/2) coordinate axis goes to the +z^,=o coordinate axis. It corresponds 
to that the point 5 moves to 6 in the RD with = 0, i.e., the point 4 in the 
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RD with = rr/2 moves to 6 in the RD with = 0. In the process, the orbit 
passes through the +1^=4^/9, +t / </>=7 7 r/i8) • • +v</>=n/i8 coordinate axes, while 
the point 4 moves in the RDs with = 4tt/9, = 7n/18, • • •, = tt/9 and 
then the point 6 moves in the RD with = tt/18. The phase of unstable 
manifold advances by an angle close to — tt. The orbit starting from the +z^=o 
(+fJ></>=n/2) coordinate axis goes into the first quadrant of configuration space 
and back to the +^=o (+^0=^/2) coordinate axis. It corresponds to that 
the point 6 moves to 1 in the RD with = 0, i.e., the point 5 moves to 
6 in the RD with = tt/2. In the process, the orbit passes through the 
+/^ =4vr/ /g coordinate axis two times, while the point 5 goes into the RD with 
= 4tt/9, then leaps to 6 and comes back the RD with = tt/2. In the 
moving processes of the point 5 from the RD with = tt/2 into the RD with 
= 4tt/9 and of the point 6 from the RD with = 4rr/9 into the RD with 
= tt/2, the directions of unstable manifolds are almost invariant. However, 
in the leaping process from the point 5 to 6 in the RD with = 4tt/9, an 
unstable direction is approximately reversed. The phase of unstable manifold 
advances by an angle close to — tt. 

So, in the periodic process, the unstable manifold returns to its original 
position, as well as the phase of unstable manifold advances by -6ir. The 
rotation number of UPO (4) is -3. 

(ii) The UPO (5) displayed in the configuration space of Fig. 4(b) passes 
through the origin. Its right limit orbit has a 13 symmetry, but not p, tt, p, 
T, <7 M , o~ u and <t 2 4 symmetries. Similarly, in the periodic process, the phase 
of unstable manifold advances by -Git. The rotation number of UPO (5) is 
-3. 

(iii) The UPO (14) displayed in the configuration space of Fig. 4(c) has 
a u symmetry, but not p, tt, p, T, <r M , o"i 3 and <7 2 4 symmetries. Similarly, in 
the periodic process, the phase of unstable manifold advances by -8tt. The 
rotation number of UPO (14) is -4. 

(iv) The UPO (15) displayed in the configuration space of Fig. 4(d) has 
p, tt, p symmetries, but not T, cr^, <t„, <7i 3 and (T 2 4 symmetries. Similarly, in 
the periodic process, the phase of unstable manifold advances by -lQn. The 
rotation number of UPO (15) is -8. 

Thus, in the rotation of unstable invariant manifolds around UPOs, we 
have determined rotation numbers. At the same time, we have also obtained 
that the advanced phase \6\ of unstable invariant manifold in a Poincare 
mapping does not exceed tt. In the Sect. V, we will present a method to 
calculate the rotation numbers in a Poincare section. 
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5 Rotation of unstable manifolds in a Poincare 
section 



According to the natural ordering in the lifted space and the occurrence of 
tangencies of manifolds, we have the region partition in the RD with symbols 
(L , R , Ri, R 2 and L 2 ) and the ordering for forward sequences[13] 

• L < »R < •R 1 < •R 2 < »L 2 . (3) 

The forward mapping preserves the ordering in regions of «Lo and »L 2 , but 
reverses the ordering in regions of »Ro, *Ri and •R 2 . 

In the RD, some symmetries of the Hamiltonian (2) are reduced, it can be 
reflected by the relation of orbit periods to sequence ones. So, we can firstly 
calculate the advanced phases of unstable directions in rotation around UPOs 
in the RD and then add the contribution of symmetries to determine rotation 
numbers. The 4 UPOs with different symmetries in the Sect. IV are still 
taken as examples. 

(i) In Fig. 5(a), we draw the periodic points encoded by R^RiR^Ri and 
stable and unstable invariant manifolds passing through the points. The pe- 
riodic points are denoted by circles. In order to illustrate the evolution of 
unstable direction around the periodic points, we take another initial point 
near the periodic point 1 on the unstable invariant manifold. The forward 
mapping of the point is also drawn in the figure and denoted by crosses. The 
arrows from periodic points to their neighboring points on the unstable in- 
variant manifolds display unstable directions. In the forward mapping from 
the periodic point 1 to 2, the symbolic sequence •R$RiR 2 Ri is shifted to 
•RqRiR^RiRq and the original unstable direction is approximately reversed. 
Since the rotation is clockwise, we obtain — 7r rotation of the unstable direc- 
tion. In the forward mapping from the periodic point 2 to 3, the symbolic 
sequence •RqR\R 2 RiRq is shifted to •RiR^RiRq and the unstable direction 
is approximately reversed. We also obtain — 7r rotation of the unstable direc- 
tion, i.e. — 2tc rotation of the original unstable direction. In the same way, — n 
rotation of the unstable direction is obtained in the forward mapping from 
the periodic point 3 to 4. In the periodic point 4, we take another neigh- 
boring point denoted by a triangle to replace the point denoted by a cross. 
In the same way, 3 x (— n) rotation of the unstable direction is obtained in 
the forward mapping from the periodic point 4 to 5, from the periodic point 
5 to 6 and from the periodic point 6 to 1. Thus, during the mapping in 
the sequence period, the original unstable direction goes back and the total 
advance of phase is — 6n. Since the orbit period is equal to the sequence one, 
i.e., the UPO (4) has not the p, ir, p symmetries, the rotation number of 
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UPO encoded by R 2 R l RlR l is -3. 

(ii) In Fig. 5(b), the periodic points with the right limit encoded by 
LqRqRiL 2 R\Ri and the stable and unstable invariant manifolds passing 
through the points are drawn. In the forward mapping from the periodic 
point 1 to 2, the symbolic sequence •LqRqRiL 2 R\Ri is shifted to •RqRi^R^RiLq 
and the original unstable direction is approximately preserved. In the for- 
ward mapping from the periodic point 5 to 6, the same result is obtained. 

In other forward mappings, 6 x (— 7r) rotation of the unstable direction is 
added. Thus, during the mapping in the sequence period, the original un- 
stable direction goes back and the total advance of phase is — 6n. Since the 
orbit period is equal to the sequence one, the rotation number of UPO (5) 
encoded by L _Rq_R 1 L 2 -R|-Ri is -3. 

(iii) In Fig. 5(c), the periodic points encoded by LqRiR 2 R\L 2 R\RqR\ and 
stable and unstable invariant manifolds passing through the points are drawn. 
In the forward mapping from the periodic point 1 to 2, the symbolic se- 
quence •LqRiR 2 R\L 2 RiRqR\ is shifted to •RiR 2 R 2 L 2 RiR R 2 L and the 
original unstable direction is approximately preserved. In the forward map- 
ping from the periodic point 6 to 7, the same result is obtained. In other 
forward mappings, 8 x (— ir) rotation of the unstable direction is added. Thus, 
during the mapping in the sequence period, the original unstable direction 
goes back and the total advance of phase is — 8n. Since the orbit period 
is equal to the sequence one, the rotation number of UPO (14) encoded by 
LqR\R 2 R\L 2 R\RqR\ is -4. 

(iv) In Fig. 5(d), the periodic points encoded by L$R\Rq and stable and 
unstable invariant manifolds passing through the points are drawn. In the for- 
ward mapping from the periodic point 1 to 2, the symbolic sequence •LqR^Rq 
is shifted to •R^RqLq and the original unstable direction is approximately 
preserved. In other forward mappings, 4 x (— 7r) rotation of the unstable 
direction is added. Thus, during the mapping in the sequence period, the 
original unstable direction goes back and the total advance of phase is —An. 
Since the orbit period is 4 times of the sequence one, i.e., the UPO (15) has 
the p, ii, p symmetries, the rotation number of the UPO encoded by LqR\Rq 
is -8. 

In the above examples describing the rotation of unstable directions around 
periodic points, the forward map corresponding to the shift with Lq or L 2 
(Ro or Ri or R 2 ) approximately preserves (reserves) the original unstable 
direction. So, we can multiply the numbers of R , R\ and R 2 in 5- letter 
symbolic sequences by one half of the ratios of orbit periods to sequence ones 
to determine rotation numbers of UPOs. 

Since the RD has the 7r-rotation symmetry, the 5-letter symbolic dy- 
namics can be reduced to the 3- letter one in the minimal domain (MD) 



9 



(s G [0,1), v G [ — 1, 1] ) [13] . The MD is partitioned and denoted by sym- 
bols L , R and Ri. The correspondence of 5-letter symbolic sequences with 
3-letter ones is L — > L , R — > R , R\ — > R 2 — > i?o and L 2 — > Lo- 
in general, the number of i?0) -Ri and i?2 in 5-letter symbolic sequences is 
twice of the number of Rq and Ri in 3-letter ones. The ratios of orbit pe- 
riods to sequence ones for the former are one half of those for the later. Of 
course, the simple repeating of 3-letter symbolic sequences in 5-letter ones 
will be removed. For example, the 3-letter symbolic sequences RqR\ and 
LqR^Ri correspond to the 5-letter ones Ro(R2)Ro{R2)RiR2{Ro)R2{Rq)Ri 
and L (L2)Ro(R2)Ro(R2)RiL2(Lo)R2(Ro)R2(Ro)Ri, respectively. Rotation 
numbers of two UPOs can be determined by calculating total numbers of 
letters R and R x in 3-letter symbolic sequences and multiplying them by 
one half of the ratios of orbit periods to sequence ones. The same rotation 
numbers of the UPOs can be also obtained by using the method for 5-letter 
symbolic sequences. Thus, using the method, we extract rotation numbers 
of 38 UPOs from symbolic sequences as given in Table I. 

6 Topological Template 

After the region partition and symbolic ordering are introduced, the families 
of stable and unstable manifolds constitute curve coordinates in the RD. 
Each stable (unstable) manifold has the same forward (backward) symbolic 
sequence. The ordering on stable (unstable) manifolds is described by that 
of forward (backward) symbolic sequences [14]. In Fig. 6, two families of sub- 
manifolds divided by the partition line «Co or mC2 have the opposite ordering. 
The ordering of stable (unstable) manifolds increases monotonically from the 
left-bottom (left-top) to right-top (right-bottom) along each unstable (stable) 
manifold. 

Along each stable manifold in zones •Lq, »Rq and »R± of Fig. 6, we con- 
strict all points in the curve to a point. The point preserves the forward 
symbolic sequence and the ordering of stable manifold. So, in the left region 
of Fig. 6, the points in three zones are reduced to three lines. Connecting the 
three lines, we obtain a belt partitioned by the symbols •Cq and •-Bo, and 
denoted by the symbols «L , »Rq and »Ri as given in the top of Fig. 7(a). 
From the left to right along the belt, the ordering of forward sequences in- 
creases monotonically. The forward mapping of the left region (s G (0, 1) 
and v G [—1,1]) in Fig. 6, i.e. the right region (s G (2,3), v G [—1,1]) in 
Fig. 3(a), can be reduced in the RD. So, in the forward mapping, the zones 
•L Q and »i?o still keep in the region (s G (0,1), v G [—1,1]) encoded by 
L « and Rq», respectively, but the zone »Ri moves in the region (s G (1,2), 
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v G [— 1, 1]) encoded by Ri*. In Fig. 6, we again partition the RD and en- 
code it by corresponding backward symbols. The two zones L « and R » are 
partitioned to five zones by the lines •Bo and «Co, as well as the zone R±» 
is partitioned to three zones by the lines »B 2 and •C*2. In the stretching and 
folding processes of forward mapping, the original three zones are mapped 
into the eight zones. We still constrict all points along stable manifolds in 
each zone. Thus, the points in eight zones are reduced to five lines in the 
left region and three lines in the right region. Connecting the eight lines, we 
obtain a belt partitioned by the symbols Lq • Bq, Co • Ri, Rq • Bq, Rq • Co, 
Ri • C 2 and R\ • B2, and encoded by the strings Lq • Ro, Lq • R±, Rq • Ri, 
Ro • R , R • L , Ri • L 2 , R\ • Ri and R\ • R\ in the bottom of Fig. 7(a). 
The ordering of the bottom belt is the same as that of the top one. So, the 
forward mapping of the region (s G (0, 1), v G [—1, 1]) can be described by a 
twisting part of topological template as given in Fig. 7(a). 

In the same way, along each stable manifold in zones *Ri, »R 2 and »L 2 of 
Fig. 6, all points in the curve are constricted to a point. A belt containing the 
point is partitioned by the symbols •B 2 and »L 2 , and encoded by the symbols 
•Rx, »R 2 and »L 2 as given in the top of Fig. 7(b). From the left to right 
along the belt, the ordering of forward sequences increases monotonically. 
The forward mapping of the right region (s G (1,2) and v G [—1,1]) in 
Fig. 6, i.e. the left region (s G (—1,0), v G [—1,1]) in Fig. 3(b), can be 
reduced in the RD. So, in the forward mapping, the zones »L 2 and »R 2 still 
keep in the region (s G (1,2), v G [—1, 1]) encoded by L 2 » and R 2 », but the 
zone »Ri moves in the region (s G (0, 1), v G [—1, 1]) encoded by Ri». We 
still constrict all points along stable manifolds in each zone and obtain the 
belt partitioned by the symbols Ri • B , R\ • Co, R 2 • C 2 , R 2 • B 2 , C 2 • Ri 
and L 2 • B 2 , and encoded by the strings R\ • R±, Ri • Ro, Ri • L , R 2 • L 2 , 
R 2 • R 2 , R 2 • Ri, L 2 • Ri and L 2 • R 2 in the bottom of Fig. 7(b). So, the 
forward mapping of the region (s G (1,2), t) G [—1, 1]) can be described by a 
twisting part of topological template as given in Fig. 7(b). 

In the two twisting parts of topological template, the belts reflect ap- 
proximately the direction of unstable manifolds. Using the twisting parts 
of topological templates, we can easily calculate rotation numbers for given 
symbolic sequences. In the same way, the 4 UPOs with different symmetries 
in Sect. II are still taken as examples. For the sequence RqRiR 2 Ri encod- 
ing the UPO(4), firstly, Rq • RoRiRjRi is obtained by shifting R x • R^RxR^. 
When an arrow is put on the «i? zone of top belt in Fig. 7(a), after a forward 
mapping, the arrow is moved on the Ro»R zone of bottom belt in Fig. 7(a). 
Since the arrow rotates clockwise to its opposite direction, we count the pro- 
cess as -1. Then, R • R\R 2 R\Rq is obtained by shifting R • RqRiR 2 Ri- 
Since an arrow on the »Rq zone of top belt in Fig. 7(a) rotates clockwise 
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to its opposite direction, we count the forward mapping as -1. Repeating 
the above process, we get the total number -6 counting the forward map- 
ping in sequence period. Since the orbit period is equal to the sequence one, 
we can thus obtain rotation number of the UPO encoded by R^R\R\R\ is 
-3. For the sequence LqRqRiL 2 R 2 R\ encoding the UPO(5), after a forward 
mapping for •LqRqRiL 2 R 2 Ri, an arrow on the top belt in Fig. 7(a) moves 
parallelly on the zone L • i? of bottom belt. We count the forward mapping 
as 0. Following the same process, we get the total number -6 counting the 
forward mapping in the sequence period. Since the orbit period is equal to 
the sequence one, we can thus obtain rotation number of the UPO encoded 
by L RlRiL 2 RlR\ is -3. Similarly, for the sequence L RiR 2 R\L 2 RiRqR\ 
{L R\Rl) encoding the UPO(14) (UPO(15)), we get the total number -8 
(-4) counting the forward mapping in the sequence period. Since the orbit 
period is equal to (4 times of) the sequence one. we can thus obtain rota- 
tion number of the UPO encoded by LqRiR 2 R\L 2 RiRqR\ (LqR\Rq) is -4 
(-8). By comparing with the former computation for the 4 UPOs from the 
definition and in a Poincare section, the same results are extracted from the 
topological template. 

By combining the two twisting parts in Figs. 7(a) (b), suspension of the 
Poincare mapping, which displays the relative position of zones in the forward 
mapping from the top belt to bottom one, is obtained. A global topological 
template of the RD is constructed by connecting the suspension with a flow 
corresponding to the Poincare mapping in Fig. 8. The template preserves 
the ordering of forward sequences encoding stable manifolds in a belt and 
the same or inverse ordering of symbolic encoding in the forward mapping of 
all parts in the belt. 

7 Conclusion and discussion 

In summary, we have presented the systematic study of the evolution of in- 
variant manifolds around unstable periodic orbits and the reduction of them 
to construct a topological template in terms of symbolic dynamics for the 
diamagnetic Kepler problem. To confirm the topological template, rotation 
numbers of invariant manifolds around unstable periodic orbits in a phase 
space, which quantify the evolution, are determined from the definition and 
connected with symbolic sequences encoding the periodic orbits. Only sym- 
bolic codes, which correspond to the forward mapping with inverse ordering, 
can contribute to the rotation of invariant manifolds. By using symbolic 
ordering, the reduced Poincare section is constricted along stable manifolds 
and a topological template, which preserves the ordering of forward sequences 
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and can be used to extract the rotation numbers, is established. The rota- 
tion numbers computed from the topological template are the same as those 
computed from the original definition. 

Since unstable periodic orbits in phase space are the skeleton of the 
chaotic system, the local evolution of manifolds near unstable periodic or- 
bits can present basic features of the global evolution of manifolds in phase 
space. One of the basic features can be quantified by the rotation number 
and embedded in the topological template. 

In the semiclassical Green's function, the phase correction is related to 
Maslov indices of the UPOs[17], which has been connected with symbolic se- 
quences of unstable periodic orbits due to boundary coding[18]. The relation 
of Maslov indices to rotation numbers of unstable periodic orbits remains to 
be determined. 
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9 FIGURE CAPTION 



Fig. 1. A typical orbit and boundary of the transformed potential for the 
diamagnetic Kepler problem at e = 0. 

Fig. 2. A image (a) and preimage (b) of the strips 1 and 2 of the funda- 
mental domain (s G [—1, 7), v G [—1, 1]) in the lifted space. 

Fig. 3. 9 lines in (a) (s G (0, 1), v G [-1, 1]) or (b) (s 6 (1,2), jj 6 [-1, 1]) 
of the reduced domain and their forward mapping in the correspondent lifted 
space. The different types of lines display the relative changes between orig- 
inal positions and their forward mappings along the s coordinate axis. 

Fig. 4. UPOs with different symmetry in the configuration space and pe- 
riodic points with unstable invariant manifolds in rotation Poincare sections: 
(a) the UPO (4); (b) the UPO (5); (c) the UPO (14); (d) the UPO (15). 

Fig. 5. Periodic points encoded by (a) R^RiR^Ri, (b) L0RIR1L2R.IR1, 
(c) LqRiR2R\L 2 R\RqR\^ (d) LqR\Rq in the 5-letter encoding and their stable 
and unstable invariant manifolds. 

Fig. 6. Stable and unstable manifolds with partition lines in the reduced 
domain. 

Fig. 7. Twisting parts of topological template describing the forward 
mapping on the regions (a) (s G [0,1), v G [—1,1]) and (b) (s G [1,2), 
v G [—1,1]). For crossing of two lines in suspension of the Poincare mapping, 
the front (back) one is denoted by solid lines (short dashes or the combination 
of solid lines and short dashes). The top and bottom belts are denoted by 
solid lines or dashes depending on their positions. The connecting lines of 
two parts of the broken belt in the forward mapping are denoted by long 
dashes. 

Fig. 8. A flow of topological template of the reduced domain. The notions 
for belts are the same as Fig. 7, except that projection of several parts in the 
bottom belt on one belt in terms of forward symbolic codings is connected 
by short dashes. The flow is denoted by short dashes. 
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Table I. Rotation numbers and symbolic sequences of 
UPOs for the diamagnetic Kepler problem at e = 0. 



No 


3-lett. Seq. & 


its Period 


5-lett. Seq. & 


its Period 


Orb. Period 


Rot. Ni 


1 


Ro 


1 


Rq 


1 


4 


-2 


2 


L Ri 


2 


LqR\L<iR\ 


4 


4 


-1 


3 


RqRi 


2 


R0R1R2R1 


4 


4 


-2 


4 


zd2 zd 
K Kl 


3 


d2 T> d2 ! > 
JXqJxiJX2^(-i 


6 


6 


-3 


5 


LqRqRi 


4 


LqRq^R\Ij2R\R\ 


8 


8 


-3 


6 


LqKiKq 


4 


LiqR\R\L2R\R\ 


8 


8 


-3 


7 


LqR\RqRi 


4 


LqRxR^Rx 


4 


4 


-1.5 


8 


R Ri 


4 


zd3 zd zd3 zd 


8 


8 


-4 


9 


h r 1 


4 


zd2 zd2 

H R 1 


4 


8 


-4 


10 


LqRqLqR\ 


5 


LqRqLqR 1 


5 


20 


-6 


11 


T E>2 d2 


5 


t td2 d2 

LqHq^ 


5 


20 


-8 


12 


LqR\LqR\ 


5 


t zd z~ zd2 z~ zd z~ zd2 
ijO-Kl^2Jli-^2-Kl-^0-Ki 


10 


10 


-3 


13 


LqRiRqRi 


5 


L1QRIR2R1 


5 


20 


-8 


14 


LqR\RqR\ 


5 


ZD ZD zd2 7~ 7D 7D zd2 

LQRIR2RI-L2R1R0RI 


10 


10 


-4 


15 


T E>2 r>2 


5 


T E>2 d2 


5 


20 


-8 


16 


LqRiRqRi 


5 


Z~ ZD2 ZD ZD T ZD2 ZD ZD 

LqK 1 KqKiL2H 1 1%2Ki 


10 


10 


-4 


17 


RqRi 


5 


RqRiR^Ri 


10 


10 


-5 


18 




5 


d3 d2 
H H 1 


5 


20 


-10 


19 


RqRiRoRi 


5 


RqR\R2R\ 


5 


20 


-10 


20 


RqRi 


5 


zd2 zd3 zd2 zd3 

RqR^R2R\ 


10 


10 


-5 


21 


RqRiRqR^ 


5 


RqRi i?2 R2R1 Ro R\ 


10 


10 


-5 


22 


LqRqLoRiRoRi 


6 


LqRqLqR\R2R\ 


6 


12 


-4 


23 


LqRqLqR\ 


6 


LqRqLqR\Ij2R2L2R\ 


12 


12 


-4 


24 


LqRqLqRi 


6 


LqRqLqR 1 


6 


6 


-2 


25 


LoRqRi 


6 


L q Rq R 1 L 2 R\ R 1 


12 


12 


-5 


26 


t d3 d2 
YjO-f«-0' tt l 


6 


T d3 d2 


6 


z? 

O 


-2.5 


27 


LqRqRI 


6 


LqR^R\Ij2R\R\ 


12 


12 


-5 


28 


L R 1 L R 1 R R 1 


6 


LqRiL2RiRqR\L2R\LqR\R2R\ 


12 


12 


-4 


29 


LqRiRq 


6 


L0RIR2L2R1R0 


12 


12 


-5 


30 


LqRiRqRx 


6 


LqRiRqRi 


6 


12 


-5 


31 


LqRiRqRiRqRi 


6 


LQR1R2R1R0R1L2R1RQR1R2R1 


12 


12 


-5 


32 


LqRiRq 


6 


LqR\Rq 


6 


6 


-2.5 


33 


LqRiRq 


6 


LqR^R^L2R\R^) 


12 


12 


-5 


34 


RqRi 


6 


RqRiR^Ri 


12 


12 


-6 


35 


RqRi 


6 


RqR\ 


6 


12 


-6 


36 


RqRi RqRi 


6 


RqRiR2Ri 


6 


12 


-6 


37 


RlR\ 


6 


RqR^ R2R1 


12 


12 


-6 


38 


RqRi 


6 


RqRi 


6 


12 


-6 
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